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Abstract 

We consider the scenario in which the knee in the cosmic ray spectrum is 
due to a change in the escape mechanism of cosmic rays from the Galaxy from 
one dominated by transverse diffusion to one dominated by drifts. We show 
that this scenario explains not only the changes in spectral slope at the knee 
and at the second knee, but can also account for the main characteristics of the 
observed energy dependent anisotropy amplitude and phase of first harmonic 
in the energy range between 10 15 and 10 18 eV. This provides a useful handle to 
distinguish this diffusion/ drift model from other scenarios proposed to explain 
the knee in the spectrum. 



1 Introduction 

A puzzling feature of the cosmic ray (CR) spectrum is the so-called knee, i.e. a 
steepening taking place at an energy E^ nee ~ 3 x 10 15 eV. Several models have been 
proposed so far in order to explain this feature, although none of them has managed 
to become broadly accepted. Some proposals focus on a possible crossover between 
different acceleration mechanisms below and above the knee [1, 2, 3], or exploit 
the possibility of a change in the particle acceleration efficiency [4, 5, 6]. Other 
scenarios include the effects of nuclear photodisintegration processes at the sources 
[7, 8], the recent explosion of a single source [9], and increased leakage from the 
Galaxy due to a change in the confinement efficiency of CRs by galactic magnetic 
fields [10, 11, 12, 13]. However, considering the information on the spectrum alone 
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it is difficult to discriminate among the predictions worked out from the proposed 
models, and so other kind of observations, namely the CR mass composition and 
the anisotropy measurements, are decisive in testing the different proposals. 

For instance, the scenario that considers nuclear photodisintegration at the 
sources predicts a CR composition that becomes lighter above the knee, while other 
(rigidity dependent) explanations predict the opposite trend. Unfortunately, the ob- 
servational status is far from being clear in this respect [14] , with some observations 
[15] suggesting that the CR composition could turn lighter while others finding that 
the heavier components become dominant [16]. 

On the other hand, anisotropy observations around and above the knee appear 
in this context as a valuable means of imposing further restrictions on the possi- 
ble scenarios. In particular, scenarios in which the knee is due to a change in the 
acceleration efficiency at the sources [4, 5, 6], or to the effects of nuclear photodis- 
integration in them [7, 8], are not expected to lead to any particular change in the 
behavior of the observed anisotropics in correlation with the observed change in the 
spectral index at the knee energy, while scenarios in which the knee is linked to a 
change in the efficiency of CR confinement in the Galaxy [10, 11, 12, 13] can predict 
such correlation. Moreover, to account for the observed energy dependence of the 
anisotropy amplitude and phase of first harmonic in the whole range of energies 
between the knee and the ankle is a highly non-trivial requirement for any of these 
scenarios. 

The aim of this work is to perform a more detailed comparison between the 
anisotropy observations and the predictions within the scenario in which the knee 
is due to a crossover from a regime in which the CR transport is dominated by 
transverse diffusion to one dominated by drift effects [12], recently reconsidered 
in detail in refs. [13, 17], and hereafter referred as the diffusion/drift model. This 
analysis will be done for the whole range of energies in which drift effects are relevant 
(E > E kn ee) and where CRs are dominated by the galactic component (E < E ank i e ~ 
5 x 10 18 eV). The CR transport in the Galaxy is determined by solving numerically 
a diffusion equation that takes account of the regular galactic magnetic field as well 
as a random turbulent field with a Kolmogorov spectrum of inhomogeneities, while 
for the extragalactic component an isotropic flux is adopted. Then, under different 
assumptions concerning the galactic field model considered, the particular set of 
parameters adopted within the given model (i.e. the field amplitudes and scale 
heights assumed for both the regular and random field components), the source 
distribution and the spectrum of the extragalactic (isotropic) flux, the results are 
contrasted to the observed CR energy spectrum, the anisotropy amplitude and phase 
of first harmonic measured by different experiments. It is found that the results show 
the right tendency to naturally account for the observed features, and hence these 
observations provide further support to this model. 



2 



2 The diffusion/drift model 



To understand the main physical ingredients of the diffusion/drift model considered 
in refs. [12, 13, 17], let us assume the CR transport in the Galaxy to be described 
by a steady-state diffusion equation V • J = Q, where Q is the source and the CR 
current is related to the CR density N through 

J = —D±V±N — D||V||iV + D A h x ViV , (1) 

with b being the unit vector in the direction of the regular magnetic field B , i.e. 
b = B /|B |, and Vy = b(b • V), while Vj_ = V — Vy. The components of the 
diffusion tensor are D\\ (along the direction of B ), D±_ in the perpendicular direction, 
while D A is associated to the antisymmetric part and determines the drift effects. 
Assuming for simplicity that the regular magnetic field is directed in the azimuthal 
direction and that in a first approach both the Galaxy and the CR sources can 
be considered to have cylindrical symmetry, one finds that D\\ plays no role in the 
diffusion equation, which involves then only and D A . These diffusion coefficients 
are obtained from the approximate expressions given in ref. [13]. Indeed, defining 
the turbulence level rj = B^ and /(B^ and + Bq), where B ran d is the rms amplitude of 
the random component of the magnetic field, the relevant diffusion coefficients are 
well described in the range of interest of this work by 

D ± ~ 8.7 x 10 27 e 3 - 24 " (1 - V ) 1/6 ( r -±) — (2) 

\.1»7 s 

and 

where r L is the CR Larmor radius. The spatial distribution of sources will be 
specified below in Section 3. 

The change in the spectrum from the diffusion to the drift dominated regimes 
can be simply understood from the fact that, as turns out from Eqs.(2) and (3), the 
energy dependence of the diffusion coefficients is given by D±(x) ~ D±(x)(E/ Eq) 1 / 3 
and D A (x) ~ D A (x)E/E [12] 1 . Hence, one finds that at low energies (E < Z E knee , 
where Z is the CR charge) the transverse diffusion is the dominant process affecting 
the transport of CRs and leading to dN/dE oc (D° ± /Dj_)dQ/dE oc E~ /3 ^ 1 / 3 , where 
(5 is the spectral index of the source, while at high energies (E ^> Z E knee ) one 
has instead dN/dE oc (D°J D A )dQ/dE oc E^- 1 . Thus, we see that for a constant 

-'Jn general, for a spectrum of random magnetic field inhomogcneities satisfying d_B 2 /dfc oc k~ a 
(where a = 5/3 corresponds to the Kolmogorov case, or a = 3/2 for a Kraichnan [18] hydromagnetic 
spectrum), one has oc E 2 ~ a . 
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source spectral index (3 ~ 2.4, one has the correct spectral slope below the knee 
(dN/dE oc E~ a , with a ~ 2.7). Each CR component of charge Z starts to be af- 
fected by drifts at an energy E ~ Z E knee , and its spectrum progressively steepens, 
with the spectral index finally changing by Aa ~ 2/3 in a decade of energy. The 
envelope of the total spectrum obtained by adding together the different nuclear 
components nicely fits the change from a spectrum oc E~ 2,7 below the knee, to one 
oc E~ 3 above it. Moreover, since all the lighter components are strongly suppressed 
above 10 17 eV, the dominant iron component will progressively steepen its spectrum 
until the overall spectrum becomes oc _E _2 - 7 ~ 2 / 3 above a few xlO 17 eV, hence also 
reproducing the behavior observed at the so-called second knee, namely a second 
steepening taking place at E s k — 4 x 10 17 eV, in which the spectral slope is close 
to a ~ 3.3. Furthermore, we have shown [13] that this scenario predicts a CR 
mass composition which is compatible with those experiments that report a heav- 
ier composition above the knee. Concerning the anisotropy observations, we have 
also shown [17] that a few % anisotropy with an excess from a direction near the 
galactic center and a deficit near the galactic anticenter direction around ~ 10 18 eV 
[19] could be accounted for by diffusion and drift currents affecting the galactic CR 
component. In this way, this simple scenario that just takes into account well es- 
tablished properties of the propagation of charged particles in regular and turbulent 
magnetic fields can provide an explanation of quite diverse galactic CR observations 
without requiring any additional assumptions. 



3 The cosmic ray spectrum and anisotropics 

The contribution to the anisotropy from galactic CRs of a given charge Z { will be 
[20] 



where Ji is the CR current corresponding to the given component, and is hence given 
by an expression analogous to Eq.(l) but involving the densities iVj corresponding 
to the CRs with charge Zj. The total anisotropy will then be given by 5 = J2i /A, 
where = Ni/N are the fractional abundances of all CR species (with N the total, 
galactic plus extragalactic, CR density). If we assume that the extragalactic flux 
is isotropic, it can be shown that this flux, which presumably gives the dominant 
contribution above the ankle, will not be enhanced by the diffusion process and hence 
will remain isotropic [17, 21]. This implies that the presence of the assumed isotropic 
extragalactic component just reduces the fractions /j of the galactic components, 
having the effect of suppressing the growth of the total anisotropy amplitude as the 
energy of the ankle is approached. 

Notice that under the assumptions of cylindrical symmetry and that the regular 
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magnetic field is in the azimuthal direction, one has that Vj_iV = ViV and also 
that the CR current will be perpendicular to the regular magnetic field (i.e. lying 
on the r — z plane). Moreover, the contribution to the CR current arising from 
the transverse diffusion will be in the direction of ViV, while the drift part will be 
orthogonal to ViV. 

As we mentioned before, at low energies the diffusion is completely determined 
by D±, and hence the contribution to the anisotropy from CRs of charge Zj is 
Si oc D^VNi/Ni, increasing with energy roughly as Si oc E 1 ^ 3 . However, at higher 
energies there is a crossover to the drift dominated regime, leading to the behavior 
Si oc E. At even higher energies, above 10 18 eV, the CR density gradually starts 
to be dominated by the extragalactic component, and hence the anisotropy will be 
determined by the intrinsic anisotropy of the extragalactic component (which for 
simplicity was assumed to be negligible in this work). 

The field models considered in this work are similar to those studied in ref. [13] 
(except for a minor modification of the vertical profile function). We will consider 
for simplicity that the regular magnetic field is directed in the azimuthal direction 
and that the system possesses azimuthal symmetry. The propagation region of the 
galactic CRs will be taken as a cylinder of radius R = 20 kpc and height 2H, with 
a halo size H = 10 kpc. Adopting cylindrical coordinates (r, 0, z) throughout, the 
regular galactic magnetic field can be expressed as Bo = (BQ tsk + BQ alo )(f), where 
B-Q lsk and B\ al ° are ^-independent functions that correspond to the disk and halo 
regular field components, respectively. The disk component is assumed to be given 

by 



«T (r, A = ^ '-^f sin (^) th> (^) , (5) 

where r c = 4 kpc is a core radius which smoothes out near the galactic center 
the overall 1/r behavior usually considered, Zd is the disk scale height, and the 
value of Bd is chosen such that BQ tsk {r b Sl z b s ) = —1 AtG at our observation point 
(fobs = 8.5 kpc, z b s = 0), the minus sign arising from the fact that the local galactic 
magnetic field is nearly directed in the —(f) direction. The th 2 (r) factor has been 
introduced just to insure that there are no singularities in the r = axis, since 
otherwise a singular drift current along that axis would be artificially produced. 
Notice also that this regular field component has reversals in its direction for r = 4, 
8, 12 and 16 kpc. 

The halo component can be described either by a field structure with radial 
reversals and also symmetric with respect to the galactic plane (denoted by R-S), 
or by a model without radial reversals (denoted NR) which may correspond to an 
independently generated halo field. In the latter case four different configurations 
arise from choosing the relative orientation of the halo field relative to the disk 
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one, and from the symmetry with respect to the galactic plane (S models having 
B(r,—z) = B(r,z), while A models having B(r,—z) = —B(r,z)). In this work, 
we will consider only the NR-A+ and NR-S+ cases, the + sign denoting that 
sign(i?g aio • BQ lsk ) > at our galactocentric radius for z > 0, since a detailed 
discussion comparing the different field configurations has already been worked out 
previously [13]. 

The halo component is then taken as 



where the normalization factor is assumed to be Bh = —1 A*G, and a halo scale 
height Zh is introduced in the vertical profile function. The function 1Z is just unity 
for the symmetric halo case, and it is taken as TZ(z) = tanh(z/z d ) in the case 
of a halo model antisymmetric with respect to the galactic plane. The Faraday 
rotation measures of pulsars and extragalactic radio sources might indicate that the 
global field structure is best described by a bi-symmetric spiral field in the disk with 
reversed direction from arm to arm, and an azimuthal field in the halo with reversed 
directions below and above the galactic plane, which could correspond to an AO 
dynamo field configuration [22]. In that case, the NR-A+ field model may reflect 
appropriately the main features of the large-scale galactic field structure. However, 
according to present observational data, the origin of the galactic field in the halo 
could also be related to the origin of the disk field, for instance if a galactic wind 
is present, somehow extending the field disk properties into the halo. Hence, the 
NR-S+ field model, as well as other variants discussed in ref. [13], are also plausible. 
For the random component, its intensity is taken as 



where z r is the corresponding scale height and B r the local rms amplitude. 

The source term in the diffusion equation is assumed to take the form Q(r, z) = 
Q(r)9(h s — \z\), with h s = 200 pc, and where the radial profile of the source distri- 
bution is considered to be either localized at a given radius r s (i.e. Q(r) = S(r — r s )), 
or as being constant within an inner and an outer radius. 

Since we deal here with a rigidity-dependent scenario, the results obtained de- 
pend on E/Z and a convolution of the contribution of each component of charge Z 
has to be made summing over all components in order to get the all-particle pre- 
diction. The source spectrum of a given CR component of charge Z is assumed to 
have an energy dependence dQz/dE oc E~^ z , where the normalization factors and 
spectral indices were obtained from the experimental data at low energies below the 






(7) 
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knee [23, 24], which provide the local observed spectral indices otz = fiz + 1/3, and 
extrapolated to higher energies without considering upper cutoffs. Also the rela- 
tive abundances of the different elements were obtained from the same fits to the 
low energy observed CR densities [23, 24]. In addition, an isotropic extragalactic 
component given by 

is added to the galactic flux as well. The energy dependence assumed for the extra- 
galactic CR contribution is similar to that inferred for the production at the sources 
(i.e. P ~ 2.4), since diffusion effects cannot enhance this flux. The amplitude of the 
extragalactic component is normalized to fit the observed fluxes for E ~ 10 19 eV. 

In figure 1, the CR spectra and anisotropy amplitudes obtained for both the 
NR-A+ and NR-S+ field models are compared to the observations from several 
experiments [25, 26, 4] 2 . The set of field parameters (B r , z d) z h , z r ) used for each 
model is indicated in the figure in /iG and kpc respectively. The effects of vary- 
ing either the field model parameters or the source distribution have already been 
discussed in detail in ref. [13]. Very briefly, a turbulence increase, which can be 
performed either by increasing the random field rms amplitude B r or by increasing 
the random field scale height z r , or also by decreasing the vertical scale heights Zd 
and Zh of the regular components, tends to suppress the drift effects, thus yielding 
larger fluxes and shifting the bent in the spectrum to larger energies. On the other 
hand, the drift effects are stronger the farther the source is located, since the drift 
currents can remove the propagating CRs from the galactic plane all the way from 
the source to the observer. Moreover, extended source distributions tend to pro- 
duce flatter spectra due to the smaller associated density gradients. The plots of 
figure 1 also show that, when considering an extended ring of sources with constant 
strength, the dominant observed flux is produced near the border which is closest 
to the observation point. 

From figure 1 we observe that, assuming quite plausible field model parameters, 
the diffusion/drift scenario is capable of reproducing remarkably well not only the 
observed CR spectrum, but also the anisotropy measurements. Notice from Eqs.(l) 
and (4) that |ViV| plays a crucial role in the enhancement of the anisotropy am- 
plitude, and thus it turns out that close by sources produce larger anisotropics as 
compared to farther ones. The effect of considering extended source distributions 

2 The anisotropy experimental data points are obtained from the amplitude of the first harmonic 
in right ascension A, which is related to the anisotropy 5 through A = S cos d cos A, where d is 
declination at which the observations are performed and A is the latitude of the direction where 
the flux is maximum. The plotted points correspond to A/ cosd for the different experiments, and 
thus are actually a lower limit for 6, as cos A cannot be recovered from right ascension harmonic 
analysis. 
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Figure 1: CR spectra and anisotropy amplitudes computed for different galactic 
magnetic field models and defining parameters, which are indicated as (B r , zj, Zh, z r ) 
in fxG and kpc respectively. The dotted straight lines correspond to the extragalactic 
flux given by Eq.(8). Also shown are the relevant experimental data points. 



8 



may be thought as given by the sum of the results produced by localized sources, 
weighted by the relative contribution of each one to the total CR density observed. 
Hence, it can happen that just a small region within the extended distribution gives 
the overwhelming contributions to the total CR flux and anisotropy observed. How- 
ever, it can also happen that different regions in the source distribution contribute 
with CR drift currents directed along very different directions and comparable in 
magnitude, giving some regions of suppressed anisotropy, while producing anisotropy 
enhancements in others. Anyhow, the results shown for localized sources should be 
indicative of the actual trend that would be obtained for different source distribu- 
tions, as follows from figure 1. 

Although some results for the NR-S+ case exhibit somewhat low anisotropies, it 
should be noticed that, even in the case of localized sources, we are actually dealing 
with rings around the galactic center, since non-central point-like sources would 
break the azimuthal symmetry assumed. Removing this simplifying assumption and 
solving the CR diffusion equation for a point-like nearby source, a further increase 
in the anisotropy amplitude might be achieved. Indeed, it has been recognized 
that the effect of discrete nearby sources, as for instance pulsars and supernova 
remnants such as Vela, Loop III and Geminga [27, 20], could actually give the 
dominant contribution to the observed anisotropies at low energies below the knee. 
The precise estimate of this effect requires however the knowledge of the spatial 
distribution, power, age, and evolution of the sources. Hence, an accurate fit of 
our results to the experimental data below the knee is not really necessary, but it 
is reassuring to observe that the predictions within this scenario tend to exhibit a 
quite acceptable agreement with observations above the knee. 

In performing the sum over all components in order to obtain the all-particle 
CR flux, we have taken into account the contribution of all nuclear species from 
hydrogen to nickel. We have actually also explored the possible contribution of 
heavier elements up to uranium, following ref. [24]. While the spectral indices of 
the elements with Z < 28 are measured at low energies (below ~ I0 14 eV), and then 
the inferred source index (3z is simply extrapolated to higher energies assuming 
that it remains constant, there is no information available on the spectral indices 
corresponding to ultra-heavy elements (i.e., with Z > 28). It was proposed in ref. 
[24] to obtain them just extrapolating the possible Z-dependence of the observed 
spectral indices of the lighter elements, using a parametrization given by — az = 
A + BZ C to fit the data. While this three-parameter expression favored a non- 
linear extrapolation (with C = 1. 51 ± 0.13), also the case of a linear extrapolation 
(i.e, C = 1) was considered and found to fit well the measured spectral indices [24]. 
Assuming the linear extrapolation of spectral indices, we found that the contribution 
of ultra-heavy elements is completely negligible within the diffusion/drift scenario. 
These ultra-heavy elements would only have a noticeable effect with the extreme 
non-linear extrapolations (for which for instance the spectral index for uranium 
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Figure 2: Anisotropy amplitudes corresponding to the main galactic CR compo- 
nents, namely protons and nuclei of helium, oxygen and iron. For comparison, also 
the total (galactic plus extragalactic) anisotropy amplitude is shown. The data cor- 
respond to the r s = 3 — 6 kpc extended source ring, the NR-A+ case and the same 
field parameters as in figure 1. 

results ajj ~ 1.9), but in our model this could not lead to a sufficient suppression 
of the spectrum above 10 17 eV, unless one adopts an upper rigidity cutoff in the 
source near E ~ Z x 10 17 eV. In ref. [24] instead, an ad hoc steepening at the 
knee of Act ~ 2 was adopted, and the resulting suppression was actually too large 
above the second knee. One could also mention that the most natural expectation 
would be to have an universal rigidity source spectrum, and hence the differences in 
the spectral indices observed at low energies most likely would reflect the effects of 
spallation processes affecting heavy nuclei. In this case no particular enhancements 
of ultra-heavy nuclei above the knee should be expected. 

The anisotropy amplitudes of different galactic CR components are plotted in 
figure 2 as functions of the energy, for the r s = 3 — 6 kpc extended source ring, the 
NR-A+ case and the same field parameters as before. As expected within the diffu- 
sion/drift scenario, the onset of the drift-dominated regime is a rigidity-dependent 
phenomenon that increases the anisotropy of all galactic CR components. It should 
be pointed out that, although light elements (particularly protons and He nuclei) 
are strongly suppressed at high energies due to leakage effects, they still contribute 
significantly to the total anisotropy. For instance, He nuclei are found to contribute 
around 30 — 40% of the total galactic CR anisotropy in the whole energy range 
between 10 14 and 10 18 eV for the case shown in figure 2, even though its abun- 
dance falls below 10% at 10 18 eV. It would certainly be of great interest to resolve 
observationally the anisotropies of separate CR components around and above the 
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Figure 3: Anisotropy phase of first harmonic in right ascension vs. energy obtained 
for the same cases as in figure 1. For comparison, experimental data measured by 
different experiments is also shown. The labels on the right indicate the directions 
of maximum CR intensity (AGC: Anti Galactic Center, NGP: North Galactic Pole, 
etc.), taking into account that 5 is contained in the r-z plane. 

knee, since that could help discriminate among some of the models proposed so far. 
Indeed, KASCADE may actually be capable of performing this kind of observations 
in the near future 3 . 

In figure 3 we display the anisotropy phase of first harmonic in right ascension vs. 
energy, again compared to the experimental observations compiled in refs. [26, 4], 
for the same cases considered in figure 1. At low energies, the transverse diffusion 
dominates and the direction of maximum intensity corresponds to the galactic center, 
as expected [17]. At higher energies the phase of the anisotropy depends actually 
on the detailed geometry of the regular field adopted, but in the models considered 
the differences are not large. However, as shown in ref. [17], for other magnetic 
field models the anisotropy maximum may point towards the galactic center near 
10 18 eV. 

As a summary, in this work we considered the diffusion/drift model to study the 
propagation of galactic CRs. In agreement with previous results [13, 17], we find that 
this scenario explains remarkably well the observed spectrum from below the knee 
up to the ankle, and in particular the changes in spectral slope at the knee and at 
the second knee. Furthermore, from the very nature of this scenario, the anisotropy 
amplitude naturally experiences a crossover from the transverse diffusion dominated 
regime to the drift dominated one, which shows up around the knee region. Indeed, 

3 We thank Ralph Engel for bringing this issue to our attention. 
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this is a particular feature that very clearly differentiates this diffusion/drift sce- 
nario from other proposals in which such a change in the energy dependence of the 
anisotropy amplitude is not to be expected. Recalling some well known statements 
claiming a kind of correlation between the spectrum steepening at the knee and the 
occurrence of an anisotropy enhancement [26], one could hence expect anisotropy 
considerations to favor this kind of scenarios against other proposals. Hence, more 
detailed and precise measurements concerning anisotropy observations will certainly 
be decisive in order to test this diffusion/drift scenario, as well as the other models 
proposed so far, and could also contribute significantly to our knowledge on the 
galactic field structure and the galactic CR source distribution. 
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